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ABSTRACT 

We use results of fits to the OPAL spectral data, obtained from non- 
strange hadronic r decays, to evaluate the difference between the vector 
and axial current correlators, YIv-a{Q 2 )- The behavior of Uv-a(Q 2 ) near 
euclidean momentum Q 2 = is used to determine the effective low-energy 
constants L°q and related to the renormalized low-energy constants 
L r w and Cg 7 in the chiral lagrangian. We also investigate how well two- 
loop chiral perturbation theory describes Uv-a{Q 2 ) as a function of Q 2 . 
This is the first determination of L^q and Cg? to employ a fully self- 
consistent model for the violations of quark-hadron duality in both the 
vector and axial channels. We also discuss the values of the coefficients 
C§y-A and Csy-A governing the dimension six and eight contributions 
to the operator product expansion representation of Hv-a{Q 2 )- 



t Permanent address: Department of Physics and Astronomy, San Francisco State University, San 
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I. INTRODUCTION 



Recently, we reanalyzed the OPAL spectral function data for non-strange 
hadronic r decays [1], the main aim being a determination of a value for the strong 
coupling at the r mass, a s {m 2 .), with a complete error analysis [2, 3]. Among 
the new elements in this analysis were the use of spectral-function moments with 
a good perturbative behavior [4], and a complete and self-consistent treatment of 
non-perturbative effects [3, 5]. This, in turn, requires a quantitative treatment of 
quark-hadron duality violations (DV) due to the clear presence of hadronic reso- 
nances in the spectral function data. The latter was accomplished by employing a 
model developed in Refs. [6, 7] that we will also use in the present article. While the 
analysis necessarily relies on this model, we demonstrated that the complete theoreti- 
cal parametrization of the spectral-function moments including the DV part provides 
a very good description of the experimental data. We chose to use OPAL data, rather 
than ALEPH data [8] because of the incompleteness of the data correlations [9] for 
the latter. 

While the central results in Refs. [2, 3] were based on fits to only the vector channel 
data, we also carried out simultaneous fits to the vector and axial channel data as 
a consistency check on our results. As a by-product, we thus have a quantitative 
theoretical description of the vector and axial spectral functions pv(f) and pa(0 from 
t = t min rs 1.3 GeV 2 to t = oo. This lets us evaluate dispersive integrals over 
Pv{t) — Pa{^) as a function of euclidean momentum Q quantitatively from the data. 
(For explicit expressions, see Eqs. (2.1) and (2.10) below.) This, in turn, allows us 
to extract certain low-energy constants (LECs) appearing in the chiral lagrangian, as 
well as some of the coefficients appearing in the operator product expansion (OPE), 
from the low and high Q 2 behavior of Uv-a{Q 2 ) , respectively. The determination of 
these LECs and OPE coefficients is the aim of the present article. As we will explain 
in detail below, we determine Hv-a{Q 2 ) by summing over experimental data up to 
t = Switch, and using our fitted spectral functions for t G [Switch, oo), where we will 
choose t sw i tc h € [tminjin 2 .] (m T is the r mass). 

This article is organized as follows. In Sec. II we give a brief overview of the 
necessary theory, including a rederivation of the Weinberg sum rules beyond the chiral 
limit tailored to our analysis. In Sec. Ill we explain our strategy for the numerical 
evaluation of Uv-a(Q 2 ) and other related functions from the OPAL data. In Sec. IV 
we present and discuss our results. We include an investigation of a fit of Uv-a{Q 2 ) to 
chiral perturbation theory (ChPT) to two-loop order. Our conclusions are contained 
in Sec. V. 
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II. OVERVIEW OF THEORY 



The LECs and OPE condensates this article aims to extract are all related to 
U V _ A {Q 2 ) defined by 1 

iWQ^jT * "(f + -ff('> , (2.1) 

with Q 2 the euclidean external momentum, and py (pa) the non-strange 1 = 1 
vector (axial) spectral functions summing the angular momentum J = 1 and J = 
contributions. Here and in what follows we take, for convenience, pa to be the axial 
spectral function without the contribution from the pion pole 

The difference py — Pa is constrained by the Weinberg sum rules [10]. It is useful 
to briefly review their derivation, beginning with the second sum rule, because of 
the subtleties involved at non-zero quark mass. Following Ref. [11], and showing 
contributions from the pion pole explicitly because it is not contained in p A (t), we 
write 



Jo 



dt w(t) (p v (t) - p A (t)) - V>{rnl) = / dz w(z) U v _ A (z) (2.2) 

2irt J\ Z \ =S0 

= ~<f dz w(z) U v ^ A (z) - -L <f dz w{z ) UyZ A (z) , 
Lki J\ z \ =sa 2m J\ z \ =so 

where w(t) is a polynomial in t, and where we split 

n v - A (z) = iJ™(z) + n™ A (z) (2.3) 



into the OPE and duality- violating (DV) parts, following Ref. [6]. The OPE part has 
the form 

n? p _5(-C 2 ) = E , (2.4) 

with, for three flavors [11, 12], 



C 2 y- A =-^ mu(n 2 )m d (n 2 ) (l - (j log ^ + c) ) + . . . , (2.5a) 

C,y- A =-\-flml + ... , (2.5b) 

where p is the renormalization scale, m Ujd (p 2 ) denote the running up and down quark 
masses and c is a numerical constant whose value is not required in what follows. In 
Eq. (2.5b), isospin symmetry has been assumed, and the Gell-Mann-Oakes-Renner 
relation has been used to express the product of the average light quark mass and 
quark condensate in terms of f n and m n . Contributions from higher-dimensional 
operators will be neglected. Next, in order to derive the second Weinberg sum rule, 



For our conventions, see Ref. [3]. 
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we choose w(t) = t. Expressing the DV part of Eq. (2.2) in terms of the DV parts of 
the vector and axial spectral functions [6], 

p^(t)- P T(t) = -lmU^_ A (t) , (2.6) 

71 

and evaluating the OPE part using Eq. (2.5), Eq. (2.2) can be rewritten as 

rso poo 

/ dt t (p v (t) - p A (t)) + dtt {p™(t) - pT(t)) (2.7) 

JO Js 

= 2 L m ir ( 1 + 3 ^ ) + 4^2 ( — ~ — ) rn u (s )m d (s )s , 

where we set p? = sq- This is the version of the second Weinberg sum rule we will 
employ. A similar derivation, choosing w(t) = 1, leads to the first Weinberg sum rule, 

pso poo 

/ dt {pyit) - PA (t)) + / dt (p^(t) - p™(t)) = 2fl , (2.8) 

Jo J s 

where we already dropped the correction coming from the OPE contributions to the 
right-hand side of Eq. (2.2), as these are numerically tiny for the sq of interest to us. 
Our conventions are such that /„■ = 92.21(14) MeV. 

The effective LECs and Cf^ are defined from the expansion of TI V -a(Q 2 ) 
around Q 2 = [13-15]: 

IV-a(Q 2 ) = SLf - 16C 8 tQ 2 + 0(Q 4 ) , (2.9) 

while the OPE condensates C%y-A and C^,v-a are defined from the high-Q 2 expan- 
sion (2.4). 

We will also use functions Ily^ involving additional polynomial weight factors 
w(x), defined by 

rivUQ 2 ) = [ dt w(t/s ) . (2.10) 

The weights we will consider are 

w k (x) = (1 - x) k , A; = 1,2. (2.11) 
Using the Weinberg sum rules Eqs. (2.7) and (2.8), one finds 

-8L* = Ih/_ A (0) = nfc^O) + ^ (2.12) 



«0 



4^(0) + 



4J? 
so 



1 _ 17 / a s (s ) \ m u (s )m d (s ) _ m% f 4 a s (s ) 



16tt 2 \ tt J / 2 2s V 3 t 



-Ycff 



yielding alternative ways to evaluate L^f. Similar equations can be derived for C| 7 . 
In these equations, we assumed that Uy_ A (Q 2 ) can be written as in Eq. (2.10), using 
the experimental spectral functions for t < s , and the approximation 



p v (t) - p A (t) « p° v (j) - p° v (t) , t>s , (2.13) 



above sq, cf. Eq. (2.7). This approximation involves the assumption that OPE contri- 
butions in principle present in the theoretical representation of py—pA are numerically 
tiny and can be safely neglected. We can test this assumption by evaluating the OPE 
corrections in Eq. (2.12), which in that equation appear as the terms depending on 
a s (so). Setting s — m r an< ^ using a s (m 2 )/7r ~ 0.1 and m u ^(m^) < 10 MeV, we 
find that the second term inside the square brackets is at most of order 10~ 5 . The 
term proportional to m%ot s (so)/so inside the square brackets is of order 4 x 10~ 4 at 
the r mass. For values of so down to 1.5 GeV 2 it will be larger, but even an order of 
magnitude will not affect our results below. 2 In fact, the contribution from the term 
2/ 2 m 2 r /sg to L\q itself is very small, about 2 x 10~ 5 . For our purposes, the dimension 
two and four OPE corrections to the approximation (2.13) turn out to be completely 
negligible, and it will be justified to drop the terms in Eq. (2.12) containing factors 
of a s (s ) in Sec. IV below. 3 

For the DV part of the vector and axial spectral functions, we will use the 
parametrization 

PvJ A (t) = e- s w*-W sin (a v/A + (5 V/A t) , (2.14) 

where ay/ a, Pv/a, 1v/a, and 5y/A are eight free DV parameters, which are fitted to 
moments of the experimental spectral functions. For a detailed discussion and history 
of this parametrization, see Refs. [6, 7, 17]. 



III. STRATEGY AND DATA 

We will evaluate T1 V _ A (Q 2 ) and n^(Q 2 ) using OPAL experimental data [1] 
for the spectral functions pv(t) and pA(t) for t < So = Switch, and approximating 
the difference pv(t) — Pa(^) by Eq. (2.13) for t > s = Switch, with values for the DV 
parameters from our previous fits to the data. We used adjusted OPAL data, updated 
to reflect current values of exclusive mode hadronic r-decay branching fractions, as 
described in Ref. [2]. We will choose t switch to be the upper end of OPAL bin N, 
obtaining 

ttW (rfl\ A+ , ,/7M /+ \ - PAW]) ( v 

^-V- A {Q ) = ^^^^WAswitch) ^ - ^ 2 (3.1) 

+ / at w(t/t swi tch> n2 • 

-^switch t + W 

Here At = 0.032 GeV 2 is the OPAL bin width and t[i] = (i - 1/2) At is the midpoint 
value of the iih bin; t sw itch = t[N] + At/2 = NAt. TIv-a(Q 2 ) is obtained by setting 
the polynomial weight w — 1. 



2 We will therefore also not worry about higher-order corrections in a s omitted from Eqs. (2.7) 
and (2.12) above, even though typically the perturbative expansions of the coefficients Cik.y-A 
converge slowly for the J = component. 

3 A less quantitative version of this argument appeared in Ref. [16]. 
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The simplest fits from which the DV parameters were obtained were fits to the 
separate vector and axial versions of Eq. (2.2) with w(t) = 1, using OPAL data to 
evaluate the moments 

/■so 

Iv/a(s ) = dt pv /A (t) (3.2) 
Jo 

through a Riemann-sum approximation like the one shown in Eq. (3.1), and varying 
so between a given s mi „ and m 2 . For w(t) = 1, all OPE contributions except the 
D = perturbative ones are negligible, and a fit to Iv/a(sq) thus yields a s and the 
DV parameters of the channel in question. 4 The value of s m i n was determined by 
requiring a good quality match between the experimental Iv/a(so) and fitted theo- 
retical representations, and stability of the fit parameters with respect to variation of 
s m i n . In this article, we will always choose t sw itch = Smm- 5 Our central results were 
obtained with the choice s m i n = 1.504 GeV 2 . 6 We have also used the more elaborate 
moments with weights 1 — (t/so) 2 and the "r kinematic weight" (1 — t/so) 2 (l + 2t/s ) 
inserted into Eq. (3.2); the perturbative part of all moments was evaluated using both 
fixed-order (FOPT) and contour-improved (CIPT) [18] perturbation theory. The non- 
trivially weighted moments also give access to the OPE coefficients C§y/A and Csy/A- 
Both pure vector and combined vector and axial channel fits were investigated. For 
a detailed account of all these fits, we refer to Refs. [2, 3]. The fit results employed 
here are always those from Ref. [2], unless otherwise noted. 

We have fully propagated all errors and correlations in the results we will report 
on below. In particular, the DV parameter values used in Eq. (3.1) are correlated with 
the data, and we have computed these correlations using the linear error propagation 
method summarized in the appendix of Ref. [3] (see, in particular, Eq. (A. 4) of that 
reference, which can be used to express the parameter-data covariances in terms of 
the data covariance matrix). 

IV. RESULTS 

We will begin with presenting the results for L\q and C 1 ^ as defined by Eq. (2.9), 
using Eq. (2.12) as well. After that, we will check the convergence of chiral pertur- 
bation theory by fitting the Q 2 dependence to the two-loop expressions for Tly-A 
calculated in Ref. [14]. Then, in Sec. IV C, we will revisit the dimension 6 and 8 OPE 
coefficients. 

A. and Cf£ 

Table 1 shows results relevant for L^. This LEC can be directly obtained from 
the second column using Eq. (2.9), or from the fourth or sixth column using Eq. (2.12). 
The DV parts of these integrals, corresponding to the second term on the right-hand 



4 a s was enforced to be equal in the two channels. 

5 We have explored taking Switch > Smim an d find that this leads to results fully consistent with 
the choice Switch = s min and no reduction in errors. 

6 This value corresponds to the upper end of OPAL bin 47. 





n y _^(o) 


n£Y A (o) 


n^(o) 


nfcr(o) 




n^r(o) 


1.408 
1.504 
1.600 


0.0522(10) 
0.0522(11) 
0.0523(11) 


-0.0039 
-0.0019 
-0.0001 


0.04019(88) 
0.04083(90) 
0.04081(90) 


-0.00042 
-0.00071 
-0.00066 


0.02738(72) 
0.02915(72) 
0.02916(72) 


0.00028 
0.00033 
0.00013 


1.504 


0.0522(11) 


-0.0019 


0.04081(91) 


-0.00072 


0.02916(72) 


0.00034 



TABLE 1: Values of Tl V -A, H-v-A' and U V-A at Q 2 = °- We always take the switch point 
between data and the duality- violating part of the spectral function at Switch = s m m (values 
for s m i n are in GeV 2 ). The superscript DV indicates the contribution from the second term 
on the right-hand side of Eq. (3.1). Duality violation parameters are from the fits of Ref. [2], 
Table 3. Results from fits using FOPT are shown above the double line, those from CIPT 
below. 

side of Eq. (3.1), are shown in the third, fifth and seventh columns. Note that the 
(absolute) errors become smaller with increasing k in Eq. (2.11), i.e., with more 
pinching at s min = t switch . We also note that the results are essentially independent 
of s min , and whether one chooses the FOPT or CIPT scheme for the evaluation of 
the truncated perturbative series. This is a consequence of the fact that the integrals 
are almost completely determined by the data part, i.e., the sum on the first line of 
Eq. (3.1), as can be seen from the always small contribution from the DV part of the 
integrals. We will henceforth use the FOPT results at s min = 1.504 GeV 2 . 
From Eq. (2.9) we find 

Lf Q = (-6.52 ± 0.14) x 10 -3 (from IL/_ A (0)) . (4.1) 

Using Eq. (2.12), one may also compute L\q from the other values shown in Table 1; 
the results are always consistent within errors. In fact, using n^f and Eq. (2.12), 
we obtain the somewhat more precise values: 

Lf Q = (-6.52 ±0.11) x 10 -3 (from n^(0)) , (4.2) 
= (-6.45 ± 0.09) x 10 -3 (from n^(0)) . 

These values are in good agreement with the value found recently in Ref. [19], 
except that our best error is twice as large. There are (at least) two reasons for this 
difference in errors, both of which point to the error in Ref. [19] being underestimated. 7 
First, Refs. [16, 19] used a DV ansatz of the functional form shown in Eq. (2.14) for 
the difference py Y — p° v , instead of using this form for each channel separately. That 
implies that Refs. [16, 19] used only four parameters to describe duality violations 
in V — A, whereas we use eight. The simplified four-parameter form assumed in 



7 For more comments on the comparison with Ref. [19], we refer to the Conclusion. 



Smin 




n'^(o) 


1.408 


-0.1356(47) 


0.0029 


1.504 


-0.1355(47) 


0.0016 


1.600 


-0.1356(47) 


0.0004 


1.504 


-0.1355(47) 


0.0016 



TABLE 2: Values of n'y_ A at Q 2 = obtained by differentiating Eq. (2.1) with respect to 
Q 2 . We always take the switch point between data and the duality- violating part of the 
spectral function at i S witch = Smin- The superscript DV indicates the contribution from the 
second term on the right-hand side of Eq. (3.1). Duality violation parameters are from the 
fits of Ref. [2], Table 3. Results from fits using FOPT are shown above the double line, 
those from CIPT below. 

Refs. [16, 19] would be valid if it happened, for some reason, that 7y = 7^ and 
Pv — Pa- Since we find very different values for 7^ and 7^ in our fits to both the 
OPAL data [2] and the ALEPH data [20], this condition is, however, not satisfied. 
The theoretical systematic error associated with the breakdown of this assumption 
is, of course, not included in the error estimates of Refs. [16, 19]. These comments 
remain relevant even if an ansatz of the form (2.14) gives a reasonable description 
of the difference pv(t) — PA(t) for large enough t: a model description of duality 
violations is only acceptable if it describes the resonance physics at higher energies 
in both the vector and axial channels individually. The second reason for our larger 
error is that Ref. [19] used the formally more precise, but in practice incomplete 
ALEPH data [9] . If ALEPH data with corrected correlation matrices were to become 
available, we anticipate that errors would be reduced relative to those obtained using 
the OPAL data for our fits as well. 

Values for the derivative of Uv-a{Q 2 ) with respect to Q 2 at Q 2 = are shown in 
Table 2. As one would expect, the results show the same robustness with respect to 
the various fits of Ref. [2] as those in Table 1. Using Eq. (3.1), we find 

C$ = (8.47 ± 0.29) x 10 -3 GeV" 2 . (4.3) 

This value again agrees with that found in Ref. [19], but our error is again about 
twice as large. Using the cubic doubly-pinched weight of Ref. [5] in Eq. (2.10) as was 
done in Ref. [19] does not lead to a smaller error in our case. The same comments 
about the reasons for our larger error as discussed above for L°q apply here as well. 

We have repeated the analysis presented here using fit values for the DV param- 
eters reported in Table 5 of Ref. [2], again taking all correlations into account. The 
results for L\q and Cf^ are virtually identical to those reported above. 
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B. Connection to chiral perturbation theory 



The LECs L\q and C$ , which are defined by the values at Q 2 = of Il v _ A (Q 2 ) 
and its derivative (cf. Eq. (2.9)), are connected to LECs in the order-p 6 chiral la- 
grangian through the relations [15] 

LV-a(O) = -8Lf (4.4a) 
= -8L r w (n) (l - 4(2/1, + vk)) + 16(2//. + fjL K )L r 9 {fj,) 

( 1 " log 4 + I log 4) - 8G ^> 0) 
lD7r 2 y m 2 3 m 2 J 

-32ml (CM-C r M-C r M) 
-32{2m 2 K + ml) (C r 6 M - C[M - C r 8 M) , 
-n' v _ A (0) = 16C$ (4.4b) 



4807T 2 m xv ^ s 



s=0 



1 / jj? 1 m 



1-log^ + llog^f )L r M, 



4 IT 2 / 2 ^ m 2 3 m 2 

Here the superscript r denotes the values of LECs renormalized at scale //, which 
below we will take to be ji = 0.77 GeV. 

The complete order-p 6 ChPT expression for Uv-a{Q 2 ) can be written as a func- 
tion of Q 2 in terms of the renormalized LECs L 9 10 and C[ 2 13 61 62 8 o 8i 87 usm g the 
results of Ref. [14]. 8 Choosing \i = 0.77 GeV, and using m\ '='139.570 MeV and 
m K = 495.65 MeV, the Q 2 dependence of U V -a(Q 2 ) in chiral perturbation theory to 
order p 6 takes the form 

U V _ A (Q 2 ) = -12.165^ - 32ml (C r 61 - C{ 2 - C r m ) (4.5) 
-32(277^ + ml) (C r 62 - C[ 3 - C r 81 ) - 16C r 87 Q 2 + R(Q 2 ; L r 9 ) , 

where R{Q 2 ', L 9 ) is a fully known non-analytic function in Q 2 coming from one- and 
two-loop contributions in ChPT, including one-loop contributions with a vertex con- 
taining Lg. Note that R(Q 2 \L 9 ) also depends on the scale //, even though we have 
not explicitly indicated any such dependence in Eq. (4.5), because we evaluated the 
numerical value of the coefficient of L r l0 at fi = 0.77 GeV. This implies that both 
R(Q 2 ; L 9 ) and all LECs appearing in this equation are to be evaluated at this value 
of /i. At Q 2 = 0, Eq. (4.5) yields Eq. (4.4), through the relation (2.9). 

If we fit Uv-a{Q 2 ) f° this order-p 6 expression, we can explore the range in Q 2 
for which order-p 6 ChPT is a valid approximation. Note that the order-p 6 expression 
is not linear in Q 2 , even though Eq. (2.9), which one obtains upon re-expanding the 



We do not quote those results here because of their length. 
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n v . A (Q 2 ) 

0.055 n — ■ — i — . — i — ■ — i — . — i — . — ^ — . — i — . — i — ■— ■ — i — . — i — ■— ■ — r 




0.030 u — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — u 
0.00 0.05 0.10 0.15 0.20 0.25 

Q 2 (GeV 2 ) 



FIG. 1: ChPT fits at order p 6 to U V -a(Q 2 )- The blue (dashed) curve includes Q 2 values 
up to 0.24 GeV 2 ; the red (continuous) curve includes Q 2 values up to 0.10 GeV 2 . 



order-p 6 ChPT expression for Q 2 < 4m 2 , is linear in Q 2 . 9 With m n and tuk fixed 
to their physical values, the data can, of course, not be used to separate the Q 2 - 
independent part of Eq. (4.5) into its individual order-p 4 and order-p 6 components 
without additional input. Such input might be obtained in the future from lattice 
studies employing a range of light quark masses. 

We have carried out a fit to Eq. (4.5) in terms of L r 10 and Cg 7 , using given values 
for all the other LECs on the right-hand side of Eq. (4.4). Specifically, we chose the 
central values L\(ji) = 0.00593 [23], 4m 2 (Q^/i) - C[ 2 (/x) - Q (X>) = -0.000067 and 
C£ 2 (/i) - C[ 3 (/i) - C^{fi) = [15] at fi = 0.77 GeV, 10 for which Eq. (4.4) becomes 

1$ = 1.521L£ o 0u = 0.77 GeV) - 0.000288 , (4.6a) 
C$ = C; 7 (n = 0.77 GeV) + 0.00328 GeV" 2 , (4.6b) 

where we also used m v = 547.853 MeV (the latter is needed for the evaluation of the 
loop contributions to the constants in these equations). 

Fits to ChPT at order p 6 are shown in Fig. 1. The blue (dashed) curve shows a 
fit with a maximum Q 2 value Q 2 max = 0.24 GeV 2 , while the red (continuous) curve 
shows a fit with Q 2 max = 0.10 GeV 2 . The fits were performed using points beginning 
at Q 2 = and spaced by 0.01 GeV 2 . These data, computed from Eq. (3.1), are 
strongly correlated, and not amenable to a standard x 2 fit, forcing us to perform a fit 
with diagonal inverse-squared-error weighting. 11 The full data correlations are then 



The threshold in the dispersive integral for Hv-a{Q 2 ) is 4m^., and not m£, since the it pole 
contribution was subtracted in defining pA(t). 

Errors on these values are only needed if one wishes to convert values for and into values 
for L r w and Cg 7 . For such an analysis, we refer to Ref. [15]. 
Thinning out the data does not help. 
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FIG. 2: Black points (crosses) show values of Cgf obtained from the ChPT fits as a function 
of Qmaxi the maximum Q 2 value used in the fit. The red point (diamond) at Q^ax = is 
the value of Eq. (4.3), for comparison. The blue points (filled circles) have been obtained 
from a fit to Eq. (4.5) with a term proportional to Q 4 added to it; see text for further 
details. 



taken into account in the quoted errors using the technique described in the appendix 
of Ref. [3]. 

Clearly, the blue curve does not provide a good fit, while the red curve does. 
We conclude that ChPT at this order gives a good match to Uv-a(Q 2 ) up to Q ~ 
300 MeV, which is about twice the pion mass. The ChPT fits are virtually linear, 
suggesting consistency with the extraction of and Cg from Eq. (2.9). From 
the ChPT fit corresponding to the red curve in Fig. 1, we obtain the values L[ = 
-4.08(9) x lO" 3 and C£ 7 = 3.97(18) x 10~ 3 GeV -2 , which correspond to 

LIq = (-6.49 ± 0.14) x 10" 3 , (4.7) 
Cg = (7.25 ±0.18) x 10~ 3 GeV" 2 . 

The value for L\q is completely consistent with Eqs. (4.1) and (4.2), but this is not 
the case for the value of Cgf, which is not consistent within errors with Eq. (4.3). 
The reason for this is that the value in Eq. (4.3) was obtained from the behavior of 
Uv-a(Q 2 ) near Q 2 = 0, while the value in Eq. (4.7) was obtained by a fit of Uv-a(Q 2 ) 
over the range < Q 2 < Q 2 max = 0.1 GeV 2 . Values for obtained by varying Q 2 max 
are shown as the black points (crosses) in Fig. 2. This figure shows that ChPT to 
order p 6 does a reasonable job in describing Uv-a(Q 2 ), but clearly order-p 8 effects, 
not included in the chiral fits, are present in the data. In contrast, the value of L\q is 
barely affected by varying Q 2 max \ it varies by less than the errors quoted in Eq. (4.2) 
over the range shown in Fig. 2. 

The presence of order-p 8 effects can be checked by redoing the ChPT fits, but now 
using Eq. (4.5) with an extra term +DQ A added. This is of course a phenomenological 



fit, because the order-p 8 structure is more complicated than just such a simple term. 
But Fig. 1 shows that the deteriorating quality of the fits with larger values of Q 2 max 
is due to some curvature showing up in Uv-a(Q 2 ) at larger Q 2 , and we expect this 
extra term to capture this curvature reasonably well. We show the results for Cff as 
a function of Q 2 max with this new term included in the fit as the blue points (filled 
circles) in Fig. 2. Indeed, the values for Cg 7 become much less sensitive to Q 2 nax) 
with values consistent with Eq. (4.3) over a much larger range. We also find that 
L^q does not change significantly as a consequence of this exercise: instead of the 
value in Eq. (4.7) we now obtain L\q = (—6.52 ±0.14) x 10~ 3 . The phenomenological 
coefficient D varies between 0.2 and 0.1 over the interval shown in the figure. 

The lesson of this exploration is that any values of U w and Cg 7 obtained from 
L^^q and C^j using (as in Ref. [15]) the order-p 6 ChPT relations of Eq. (4.4) must be 
treated with some care. While terms beyond order-p 6 in the chiral counting associated 
with higher powers of Q 2 can be removed by taking Q 2 to zero, those associated with 
higher powers of the quark masses are fixed by the non-zero, physical meson masses 
and cannot be removed. Such contributions are present in the relations between 
and L r 1{) and C$j and Cg 7 to arbitrarily high chiral order. The variation in the fitted 
value of Cg 7 with Q 2 max (the source of the variation of C§ 7 displayed in Fig. 2, cf. 
Eq. (4.6b)) indicates non-trivial Q 2 -dependent contributions of order-p 8 and beyond, 
raising the possibility of analogous mass-dependent, Q 2 -independent order-p 8 (and 
beyond) contributions as well. 

The impact of such order-p 8 (and higher) contributions will be more significant 
for the relation between Cg 7 and Cg 7 than for that between L\q and L r w since in 
the former case the missing order-p 8 terms are only one chiral order higher than the 
LEC of interest, Cg 7 , whereas in the latter case the missing terms begin two chiral 
orders higher than L\ . Even so, the order-p 8 and higher contributions need not be 
completely negligible for L r 1Q . In fact, ignoring the contributions of the order-p 6 LECs 
C12, 13,61,62,80,81 to the relation between L°q and L{ , the effects of the mass- dependent 
order-p 6 terms are significant, changing the coefficient of L[ from 1 at order-p 4 to 
1.521 at order-p 6 in Eq. (4.6a), and altering the best fit results for U w by about 30% 
between order-p 4 and order-p 6 . A further shift in L r w by about 0.3 2 ~ 10% due to 
order-p 8 effects would thus not be unexpected. Similarly, a difference of about 30% 
between the order-p 6 and order-p 8 values for Cg 7 would not be surprising. 

In conclusion, if estimates for L r 10 and Cg 7 obtained from and Cg 7 are used 
in the computation of some other physical quantity, propagating the error quoted 
in Eqs. (4.2) and (4.3) would not include additional systematic errors due to the 
omission of order-p 8 terms in ChPT. 

C. V — A condensates 

In this subsection, we consider the values of the OPE coefficients C 6 y~A and 
Cs,v-A, defined in Eq. (2.4). In Ref. [2] we presented fit results for C 6 y/A and 
Cgy/A obtained using sum rules involving weights up to degree three, from which it 
is straightforward to obtain C 6 y_A and C 8 y-A- From the fits at s min = 1.504 GeV 2 , 
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and including all correlations, we find the values 

C 6 y- A = (-10.5 ±2.8) x 10~ 3 GeV 6 (FOPT) , (4.8a) 

= (-11.3 ± 2.4) x 10~ 3 GeV 6 (CIPT) , 
Csy-A = (14 ± 7) x 10~ 3 GeV 8 (FOPT) , (4.8b) 

= (16 ± 6) x 10~ 3 GeV 8 (CIPT) . 

Changes as a function of varying s m i n are small compared to the errors shown in 
Eq. (4.8). 

It is interesting to compare these values with those we would obtain from the 
original OPAL data, to which no correction reflecting modern values for the r hadronic 
branching fractions have been applied. In this case, we find, using the fits reported 
in Table 5 of Ref. [3]: 

C 6tV -A = (-3 ± 4) x 10~ 3 GeV 6 (FOPT) , (4.9a) 

= (-4 ± 4) x 10~ 3 GeV 6 (CIPT) , 

C 8 y-A = (-3 ± 12) x 10~ 3 GeV 8 (FOPT) , (4.9b) 

= (0 ± 12) x 10~ 3 GeV 8 (CIPT) . 

The results for C e y_A an d Csy-A are barely consistent between the updated and 
original OPAL data. The relatively large differences between the "updated" and 
"original" data are not a big surprise: these OPE coefficients parametrize the most 
subleading part of the fits carried out in Refs. [2, 3]. Moreover, it was found that 
the fits reported in Table 5 of Ref. [2], while consistent with simpler fits, are at the 
"statistical edge" of what can be extracted from the OPAL data. 

One can avoid using the fits of Table 5 of Ref. [2] by employing the sum rule (2.2) 
with a judicious choice of the weights w(t). As we have seen, pv{t) — Pa(0 can be 
obtained from the simpler fits reported in Table 3 of Ref. [2] . An obvious possibility is 
to choose w(t) = t 2 or w(t) = t 3 , for which the right-hand side of Eq. (2.2) immediately 
yields C 6 y_A, respectively, — C 8 y-A- We find results consistent with those reported 
in Eq. (4.8), with comparable errors. 

However, using the moments of Ref. [5], which involve a double-pinching factor 
(t — ^switch) 2 , we can do better. 12 Choosing w(t) — (t — Switch) 2 or w(t) — (t — 
£ S witch) 2 (£ + 2t swit ch), Eq. (2.2) implies 

N 

Cey-A = J2 At " tswitch ) 2 (M*[*D " PA(t\i})) - 2/ 2 (m 2 r - t switch ) 2 



i=i 

"OO 



/oo 
dt (t - t switch ) 2 (p^(t) - P 7(t)) , (4.10a) 
- switch 

N 

Csy-A = ~J2 At "Switch) 2 (t[i] +2t switch ) (p v (t\i}) - PA(t[i])) 



i=i 

+2/ 2 (m 2 - t S witch) 2 (m 2 + 2t switc h) 

POD 

/ dt(t- t switch ) 2 (t + 2t switch ) (p^(t) - p° v (t)) . (4.10b) 

J ^switch 



12 This method was also employed in Ref. [19]. 
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TABLE 3: C 6 ,v-a (in GeV 6 ) and C 8 ,v-A (in GeV 8 ) from Eq. (4.10). The superscript DV 
indicates the part coming from the DV integrals in Eq. (4.10). Duality violation parameters 
are from the fits of Ref. [2], Table 3. Results from fits using FOPT are shown above the 
double line, those from CIPT below. 



In these expressions, the sums over bins, as well as the pion-pole terms, are obtained 
from data, the latter with negligible errors. 13 These sum rules have two advantages: 
(1) they suppress the data at higher t, which have larger errors, and (2) they suppress 
the contribution from the DV-integral terms [5, 19], replacing these contributions, in 
effect, by the pion-pole terms, which are known with great precision. 

We present the results in Table 3. The DV parts are significantly smaller than 
those we would obtain with w(t) = t 2 or w(t) = t 3 , especially for C^y-A, but also for 
Csy-A- And indeed, errors are also significantly smaller than those of Eq. (4.8), as 
we expected. We note, however, that there is some discrepancy between the values of 
Table 3 and Eq. (4.8). The results of Table 3 are based on results from simpler and 
more stable fits reported in Table 3 of Ref. [2]. 14 Therefore, we take as our central 
results for C^y-A and C$,v-a the values from Table 3 above, 

C 6 , y _A = (-6.6 ± 1.1) x 1(T 3 GeV 6 , (4.11) 
C 8 ,v-a = (5 ± 5) x 1CT 3 GeV 8 , 

where the central values are the averages of the values in Table 3, and the errors have 
been obtained by adding the fitting error at s m i n = 1.504 GeV 2 and the variation as 
a function of s m i n in quadrature. 

In Fig. 3 we compare our results with other results in the literature. Updating 
the OPAL data without including DVs in the analysis causes the central values of the 
OPAL-based results of Ref. [22] to shift from C %y _ A = -5.4 x 1CT 3 GeV 6 , C$ y „ A = 
-1.4 x 1CT 3 GeV 8 to C 6 ,y_A = -5.0 x 10~ 3 GeV 6 and C 8 y- A = -3.4 x 10~ 3 GeV 8 . A 
comparison of the latter set to the results of the present analysis then shows directly 
the impact of the inclusion of DVs. 15 



Ordcr-a s corrections from dimension two and four terms in the OPE are again completely negli- 
gible. 

For an extensive discussion of the quality of these fits, we refer to Ref. [2]. 

The reader should note that, for the s employed in the fits of Ref. [22], integrated DVs have the 
opposite sign to those shown in Table 3. 
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FIG. 3: Comparison with other recent values [19, 21, 22] for 6*6, v-A (left panel) and Cg,v-A 
(right panel). 



We note in particular that our values do not agree with those found in Ref. [19]. 
While our discussion above indicates that the determination of C§y_A an d C S y-A 
is limited by the quality of the data, we also recall that Ref. [19] used a much more 
restricted parametrization of duality violations in the V —A channel, with four instead 
of eight parameters. 

It is interesting to compare our results for C§y_A with an analytical expression 
that is available at the next-to- leading order [24] (see also ref. [25]): 

32 / 119 \ 
Ce,v-A = - Y 7T { 1+ 24^ as ( s °)) as ( s °^ + Ps)(gg(so)) 2 (4.12) 

2 

- 3 afOoXPi + P5)(qq(s )) 2 ■ 

The parameters pi^ and p\$ parametrize deviations from the so-called vacuum sat- 
uration approximation (VSA), in which they are all normalized to unity. Values for 
pi t 5 from our fits have already been discussed in Refs. [2, 3]. Numerically, at so ~ m 2 . 
the second line of Eq. (4.12) only contributes about a few percent, so that precise 
values for p\^ are irrelevant. On the other hand, in the VSA the first line of Eq. (4.12) 
yields 

C^ A = -4.4 x KT 3 GeV 6 , (4.13) 

where (qq(m 2 T )} = -(272 MeV) 3 [26], together with our result for a s (m 2 T ) has been em- 
ployed. As the next-to-leading order correction in Eq. (4.12) amounts to about 50%, 
an error of that size should be attributed to the numerical value (4.13). Therefore, the 
difference between our central fit result of Eqs. (4.11) and (4.12) could either be due 
to higher-order QCD corrections or a breaking of the VSA. At any rate, no significant 
deviations from the VSA are observed and the results in Eqs. (4.11) and (4.13) are 
nicely compatible. 
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V. CONCLUSION 



We used results of earlier fits to the non-strange vector- and axial-channel spectral 
functions obtained from OPAL hadronic r decay data in order to estimate the low- 
energy constant combinations and , as well as the operator product coefficients 
Cey-A and Cgy-A- Our best values are 

(-6.45 ± 0.09) x 10^ 3 , (5.1) 
(8.47 ± 0.29) x 10~ 3 GeV -2 , 
(-6.6 ±1.1) x 10~ 3 GeV 6 , 
(5 ± 5) x 10~ 3 GeV 8 . 

For a comparison with the values of L\q and Cff obtained in Ref. [19], we refer to 
Sec. IV A. For comparisons with other values for C§y-A and Csy-A obtained in the 
literature, see Fig. 3. As emphasized in Sec. IV C, for C^y-A an d Csy-A the results 
are rather sensitive to small variations in the data, and to the details of the fits. In 
contrast, we expect the results for L\q and C$ to be rather robust, since these values 
are dominated by the low-Q 2 range of the data, where the experimental errors are 
small. For a comparison of the low-Q 2 behavior of Uv-a{Q 2 ) with ChPT to order p 6 , 
we refer to Sec. IV B. We find that order-p 8 effects, not included in our chiral fits, are 
clearly visible in Cg 7 , but not in L r lQ . This is consistent with what one would expect: 
taking into account order-p 6 terms stabilizes the values of the LECs at lower order. 

We demonstrated in both Ref. [3] and Ref. [2] , that our fits satisfy both Weinberg 
sum rules, as well as the DGMLY sum rule for the pion electromagnetic self-enegy 
[27] within errors, though none of these were enforced in the fits. The situation 
is thus very much analogous to that of the analysis of Refs. [16] and [19]. There, 
the set of "acceptable" DV parameter combinations was generated by requiring the 
corresponding DV contributions to the Weinberg and DGMLY sum rules to be such 
that all three sum rules were satisfied within the experimental errors on the data 
part of these sum rules, i.e., the integral from to So in Eq. (2.2). On this point, 
there is thus no relevant difference between the strategies employed in Refs. [2, 3] and 
Refs. [16, 19]. 

There are important differences, however. First, Refs. [16, 19] started from an 
ansatz of the form (2.14) for the DV part of pv — Pa involving only four parameters 
rather than four for each of the two channels separately. The possibility that the 
vector and axial DV contributions are such as to allow the V — A combination to be 
expressed in this simplified form, however, is not supported by the results of our fits 
to the individual vector and axial channels. Furthermore, the procedure of Ref. [16], 
described in more detail in Ref. [28], does not take into account the correlations 
between the data and DV parameters induced by the use of the Weinberg and DGMLY 
sum rules. Neither were the correlations between the data and the DV parameters 
taken into account when using their results to evaluate the quantities of interest, L c ^ , 
etc. In our analysis, we have taken these correlations fully into account, and find them 
to have a significant effect. 

Finally, most of the earlier results shown in Fig. 3 are based on ALEPH data 
[8, 29]. At least for those earlier works which employed the 2005/2008 version of 
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these data [8], the incompleteness of the 2005/2008 correlation matrices [9] should 
be born in mind when appraising these results. We wish to reiterate the expectation 
that inclusive spectral functions extracted from BaBar or Belle would be of great help 
in reducing the uncertainties on C§y-\ and Cg,v-A, for the reasons already discussed 
in Ref. [2]. 
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